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Abstract. The Nagel-Schreckenberg traffic flow model shows a transi- 
tion from a free flow regime to a jammed regime for increasing car density. 
The measurement of the dynamical structure factor offers the chance to 
observe the evolution of jams without the necessity to define a car to be 
jammed or not. Above the jamming transition the dynamical structure 
factor S(k,u) exhibits for a given k- value two maxima corresponding to 
the separation of the system into the free flow phase and jammed phase. 
Analyzing the fc-dependence of these maxima the backward velocity of 
the jams is measured. We find that the jam velocity neither depends on 
the global density of the cars nor the maximal velocity of the model. 

1 Introduction 

In 1992 Nagel and Schreckenberg [0 introduced a simple cellular automata 
model, which simulates single-lane one-way traffic, in order to study the tran- 
sition from free flow traffic to jammed traffic with increasing car density. The 
behavior of the model is determined by three parameters, the maximal velocity 
ii max , the noise parameter P and the global density of cars p — N/L, where N 
denotes the total number of the cars and L the system size, respectively. The 
variables describing a car i at time t are its position Ti G {1, 2, L}, its velocity 
Vi G {0, 1, u m ax} and the gap <?i, which is the number of unoccupied cells in 
front of the car. Using periodic boundary conditions, the following update steps 
are applied in parallel for each car: 

Vi i — > min(u max ,u 1 + 1), 
Vi i — ► mm(vi,gi), 

Vi i — ► max(0, Vj, — 1) with probability P, 

n i — > ri + Vi . 

Increasing the global density for fixed P and v max jams occur above a cer- 
tain critical value p c Q . Contrary to the forward movement of all particles the 
jammed region is characterized by a backward moving of shock waves, i.e. jams 
are nothing else than backward moving density fluctuations. This property of 
jams was found in real traffic flow already in the 50 's 0. Investigations of the 
Nagel-Schreckenberg traffic flow model show that crossing the critical point a 
transition takes place from a homogeneous regime (free flow phase) to an inho- 
mogeneous regime which is characterized by a coexistence of two phases (free 
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flow and jammed traffic). Thereby, the free flow is characterized by a low local 
density and the jammed phase by a high local density, respectively. Due to the 
particle conservation of the model the transition is realized by the system sep- 
arating into a low density region and a high density region ||. Therefore, the 
dynamical structure factor is an appropriate tool to investigate the decomposi- 
tion of the two phases above the transition point. Additionally the dynamical 
structure factor has the decisive advantage that both phases could be distin- 
guished in a natural way by the sign of their characteristic velocities. Thus no 
artificial definitions of jams are needed which were used in previous investigations 
and lead to controversial results (see for instance H[|,||]g|] ) • 



2 Simulation and Results 

The dynamical structure factor is defined as follows: Let 

{1 if cell r is occupied at time t . . 

otherwise. 1 ' 

The evolution of r\ r ^t leads directly to the space-time diagram where the propa- 
gation of the particles can be visualized. Figure |l| shows a space-time diagram 
of the system above the critical value. Traffic jams are characterized by the 
backward movement of the high density regions. The dynamical structure factor 
S(k, u>) is given by 



S(k,w) = 





space -> 



Fig. 1. Space-time plot for u max = 5, P — 0.5, and p > p c . Note the separation 
of the system in high and low density regions. 
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Fig. 2. The dynamical structure factor S(k,uj) below (upper figure) and above 
(lower figure) the critical value. The ridges with maximal S^fcjtj) indicates the 
various modes. 



where the Fourier transform is taken over a finite rectangle of the space-time 
diagram of size I x T (see Fig. |lj). 

In Fig. H we plot the dynamical structure factor below and above the tran- 
sition, respectively. Below the transition S(k,u>) exhibits one mode formed by 
the ridges. The values of the modes are obtained from a determination of the 
position of the maxima of S(k, uj) for fixed ui or fc, respectively. This mode is 
characterized by a positive slope v = ^ corresponding to the positive velocity 
of the particles in the free flow phase. Increasing the global density of the cars, 
a second mode appears at the transition to the coexistence regime. This second 
mode exhibits a negative slope indicating that it corresponds to the backward 
moving density fluctuations in the jammed phase. 
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Fig. 3. The maximal values of the dynamical structure factor for t> max = 5, 
P = 0.5, and p > p c . The lines with positive slope corresponds to the free flow 
phase and the negative slope corresponds to the jammed phase, respectively. 

In Fig. H we plot the position of the ridges of S(k, ui) in the jammed regime. 
Both modes are characterized by a linear dispersion relation = const), 
i.e. each phase is characterized by one velocity only. Since the modes of both 
phases can clearly be distinguished (note the logarithmic scale in Fig. ||) and 
no other mode occur we think that this indicates that above the transition a 
real phase separation of the system takes place. This is in contrast to results of 
previous investigations which examine the steady state correlation function and 
concluded that there is no real phase separation of the system . 

A linear regression of the modes of Fig. || yields the velocities of both phases. 
In Fig. |i| we plot the velocities of the free flow and jammed phase for different 
values of P, w ma x, and p. Independent of the value of the global density, the 
free flow phase is characterized by the velocity V{ = v max — P, i.e. Vf equals 
the average velocity of the particles in the limit p — > and the cars can be 
considered as independent particles. This agrees with investigations of the local 
density distribution and the steady state structure factor (a; = 0) ||. The plotted 
values of the jam velocity uj are obtained from simulations with P = 0.5 and 
various values of p and v max . The velocity of the jams neither depends on the 
global density nor on the maximum velocity. The jam velocity is a function of 
the noise parameter only. Using the structure of S(k, lo) we are able to determine 
this P dependence of Vj for the first time and the results are plotted in Fig. H. 



3 Conclusions 

We studied numerically the Nagel-Schreckenberg traffic flow model. The investi- 
gation of the dynamical structure factor allowed us to examine the transition of 
the system from a free flow regime to a jammed regime. Above the transition the 
dynamical structure factor exhibits two modes corresponding to the coexisting 
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free flow and jammed phase. Due to the sign of their characteristic velocities Vf 
and Vj, both phases can clearly be distinguished. We investigated the dependence 
of Vf and Vj on the system parameters systematically and found that the jammed 
velocity depends on the noise parameter only The characteristic velocity of the 
free flow phase turned out to be equal to the velocity of free flowing cars in the 
low density limit (p — > 0), i.e. cars of the free flow phase behave as independent 
particles. 
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Fig. 4. The characteristic velocities of the free flow phase vf (left) and jammed 
phase Vj (right) for various values of the parameters P, v maK} and p. The dashed 
lines in the left figure correspond to the average velocity v max — P of a single 
particle. The values of Vj are obtained from simulations with P — 0.5 and various 
values of p and w max . 
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Fig. 5. The jam velocity Vj as a function of the noise parameter P. 



